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ABSTRACT 



Context. Water is a key molecule in the star formation process, but its spatial distribution in star-forming regions is not well known. 
Aims. We study the distribution of dust continuum and H 2 and 13 CO line emission in DR21, a luminous star-forming region with a powerful 
_^ outflow and a compact H n region. 
Cln 1 Methods. Herschel-HWl spectra near 1100 GHz show narrow 13 CO 10-9 emission and H 2 In - Ooo absorption from the dense core and broad 
emission from the outflow in both lines. The H2O line also shows absorption by a foreground cloud known from ground-based observations of 
low-7 CO lines. 

Results. The dust continuum emission is extended over 36" FWHM, while the 13 CO and HtO lines are confined to »24" or less. The foreground 
^2 ' absorption appears to peak further North than the other components. Radiative transfer models indicate very low abundances of ~2xl0~ 10 for H 2 
1 1 and ~8xl0~ 7 for l3 CO in the dense core, and higher H2O abundances of ~4xl0~ 9 in the foreground cloud and ~7xl0~ 7 in the outflow. 

Conclusions. The high H 2 abundance in the warm outflow is probably due to the evaporation of water-rich icy grain mantles, while the H 2 
y—{ abundance is kept down by freeze-out in the dense core and by photodissociation in the foreground cloud. 
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only probe the small traction or the gas at I 
(IVan der Tak et al.ll2006l;IWatson et al.lr2007l) . 



*-0 The water molecule is a key species throughout the formation 
O f stars and planets. In the gas phase, it acts as a coolant of col- This P a P er P resents observations of an H 2 ground state 
-> lapsing interstellar clouds; in the solid state, it acts as glue for llne at >3x hl S her an S ular resolution than previously possi- 
dust grains in protoplanetary disks to make planetesimals; and ble for such llnes - Through radiative transfer models, we com- 
> ! as a liquid, it acts as transporter bringing molecules together on P 3 " 5 the abundance distribution of H 2 with that of CO 
\Tl ' planetary surfaces, a key step towards biogenic activity. The first and dust The source DR21 < Mam ) ls a high-mass protostel- 
/S ; role is especially important for high-mass star formation which lar ob J ect (£=45,000L B ) located in the Cygnus X region at 
& , depends on the balance between the collapse of a massive gas <*=l-7kpc (Schneider et al. 2006), about 3' South of the well- 
cloud and its fragmentation (Zinnecker& Yorke 2007). " known DR21(OH) object (also known as W75S). Maps of 
Interstellar H 2 is well known from ground-based ob- the L2mm dust emission show a dense core with a mass of 
servations of the 22 GHz maser line. Previous space- 600-1000 M and a size of 0.19x0. 14pc FWHM, surrounded 
based submm and far-IR observations have measured H 2 by an extended envelope with mass 4750 Mp and size 0.3 pc 
abundances rangin g from IP' 8 in cold g as to 10 4 in dMotte etal. || 2007|) . Gas densities of 10 5 -10 6 cm j are derived 
warm gas (ISO: 'Van Dishoeck & Helmich 1996; SWAS: &om both the mm-wave continuum and HCN and HCO line 
Melnick & Bergin 2005; Odin: Bjerkeli et al. 2009) but did not emission (Kirby 2009). Signs of active high-mass star forma- 
have sufficient angular resolution to determine the spatial dis- tlon are the bn 8 ht mld " IR emission (272 Jy at 21 ^m), the pre s- 
tribution of H 2 In contrast, space-based mid-IR and ground- ence of an H 2° 22 GHz maser < see catalo S of iBraz & Epchtein| 
based mm-wave observations have high angular resolution but 1983) and emission from ionized gas extending over 20-30 

(Roelfsemaet al. 1989). Together with the powerful molecular 

* Herschel is an ESA space observatory with science instruments outflow (Garden et al.|| 199 1|) these signs indicate that the source 

provided by European-led Principal Investigator consortia and with im- is relatively evolved within the embedded phase of high-mass 

portant participation from NASA star formation, beyond the 'ultracompact H 11 region' phase. 
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Fig. 1. Spectra of 13 CO 10-9 (left) and H 2 In - o (right) lines toward DR21, taken with the WBS backend. 



2. Observations 

The DR21 region was observed with the Heterodyne Instrument 
for the Far-Infrared (HIFI; De Graauw et al, this volume) on- 
board ESA's Herschel Space Observatory (Pilbratt et al, this 
volume) on June 22, 2009. Spectra were taken in double side- 
band mode using receiver band 4b, with vlo= 1107.990 GHz and 
vrp=6GHz. The data were taken during the performance ver- 
ification (PV) phase using the double beam switch observing 
mode with a throw of 2'5 to the SW. The position observed is 
R.A. 20 :39:02.38, Dec +42:19:3 3.5 (J2000), close to radio peak 
C from iRoelfsema et al.l dl989t) . A strip map was made in the 
N-S direction, spanning offsets from +90" to -90" at a 1075 
spacing, half the beam size of 21" FWHM at our observing fre- 
quency, which corresponds to 0.17pc at the distance of DR21. 
This beam size was measured before launch and is 10% larger 
than the diffraction limit due to spillover effects. 

Data were taken with two backends: the acousto-optical 
Wide-Band Spectrometer (WBS) which covers 1 140 MHz band- 
width at 1.1 MHz (0.30kms _1 ) resolution, and the correlator- 
based High-Resolution Spectrometer (HRS), which covers 
230MHz bandwidth at 0.48 MHz (0.13kms _1 ) resolution. Two 
polarizations are available except for the HRS data of H2O. 

The system temperature of our data is 340-360 K DSB 
and the integration time is 67 seconds per position (ON+OFF). 
Calibration of the raw data onto FJ scale was performed by 
the in-orbit system (Roelfsema et al, in prep); conversion to 
T m b was done assuming a beam efficiency of 0.67 as estimated 
by the Ruze formula and validated by raster maps of Saturn 
(M. Olberg, priv. comm.). Currently, the flux scale is accu- 
rate to xT0% which will improve when the telescope efficiency 
and sideband ratio are measured on Mars. The calibration of 
the data was performed in the Herschel Interactive Processing 
Environment (HIPE; Ott 2010) version 2.1; further analysis was 
done within the CLASS3 package. After inspection, data from 
the two polarizations were averaged together to obtain rms noise 
levels of 97 mK on 0.5 MHz channels for the WBS data, 195 mK 
on 0.24 MHz channels for the 13 CO HRS data and 244 mK on 
0.24MHz channels for the H 2 Q HRS data. 



3. Results 

Figure Q] shows the WBS spectra from both receiver sidebands 
at each of the 18 offset positions. The spectra from both back- 
ends show a continuum signal and two line features. The first 
line feature is seen in emission at v — 1101.357 GHz LSB or 
1 1 14.623 GHz USB, which we identify with the 13 CO /=10->9 
line at 1 101.349597 GH2Q. This line has an upper level en- 
ergy (E„) of 291 K and a cr i tical d ensity of lxl0 6 cirT 3 , using 
collision data from iFlowerl d200ll) . The second feature has a 
mixed emission-absorption profile and lies at v = 1 102.423 GHz 
LSB or 1113.557GHz USB, which we identify with the H 2 
-/k p k„ = 1h -> Ooo line at 1113.34306GHz. This line has E u = 
53 K and a critical dens ity of 3-6xl0 8 crrT 3 , using collision data 
from lFaure et al.l d2007l) . No other lines are detected. 

Figure [2] shows the HRS spectra of the 13 CO and H 2 
lines at the central position, with the continuum subtracted. The 
13 CO line profile has a double structure which is seen to be 
well reproduced by the sum of two Gaussians: a narrow com- 
ponent (AV=4.9kms~ 1 ) centered at VLSR=-3.2kms~ 1 , which 
we attribute to the molecular cloud core (protostellar envelope), 
and a broad component (AV=15.2kms _I ) centered at Vlsr=- 
2.0kms~ 1 , which we ascribe to the molecular outflow. In con- 
trast, the H 2 line profile shows two absorption features, with 
strong emission in between, and weaker emission on the far blue- 
and redshifted sides. Overplotted is a decomposition with three 
Gaussian components: broad (AV=23.8kms~ I ) emission cen- 
tered at Vlsr— -l-0kms -1 , and absorptions at VLSR=-2.6kms~ 1 
and +7.6 km s 1 with A V=5.0 and S.Okms" 1 . This sh ape re- 
sembles the prediction by Poelman & van der Tak) (120071) for the 
H 2 ground state lines towards high-mass protostellar envelopes 
for the case of a constant, low H 2 abundance (H 2 0/H 2 ~10~ 9 ) 
with two modifications. First is the broad emission also seen in 
I3 CO which is likely due to the outflow. Second is the absorp- 
tion at VLSR=7.6kms~ 1 which is known from gr ound-based ob- 
servations of low-/ CO lines dJakob et al.l l2007) and likely due 
to a foreground cloud. This absorption is not seen in the 13 CO 
10-9 line, nor in ground-based observations of mid-7 HCN and 



1 http://www.iram.fr/IRAMFR/GILDAS 



2 S pectroscopic data are taken from the CDMS catalog (Mttller et al. 
I2005h at |http://cdms.de| 
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Fig. 2. Spectra of 13 CO 10 - 9 (top) and H 2 1 , i - 00 (bottom) 
lines toward the central position, taken with the HRS backend, 
with Gaussian decompositions overplotted. 



HCO + lines, which indicates a low temperature (~10 K) and den- 
sity ( ; 10 4 cm 3 ) for the foreground cloud. Both absorptions are 
also seen in the SWAS spe ctra of the 0-H2O a nd C° ground-state 
lines, but not in 13 CO 5^1 (lAshbv et alJl2000h . 

Figure[3]shows the spatial brightness distribution of the con- 
tinuum and the lines. The points are the results of Gaussian pro- 
file fits for the lines, and of linear baseline fits for the contin- 
uum, at each offset position. For the H2O line, separate fits were 
made for the emission and the two absorption components. The 
curves in Figure [3] are Gaussian fits to the observed spatial dis- 
tribution. The dust emission peaks at an offset of (8.7+0.5)" 
North of the nominal position, and the emission distribution 
has an FWHM width of (35.7+0.5)", consistent with ground- 
based measureme nts at longer wavelengths (Gibb etaTI 120051; 
iMotte et all 120071) . The peak of the 13 CO emission is at off- 
set (9.65+0.5)", consistent with the dust peak within the com- 
bined error, but the FWHM of the 13 CO emission is (24.1+0.4)", 
essentially unresolved. The FbO emission component has an 
FWHM of (25.2+2.5)", similar to 13 CO, but its peak is shifted 
further South, at offset (6.4+1.1)". The spatial distribution of 
the H2O absorption at V - -3kms _1 is the same as that of the 
H2O emission, peaking at offset (6.7+1.6)" with an FWHM of 
(26.2±3.3)". In contrast, the H2O absorption at V - +7.5kms~ 1 
peaks much further North, at offset (14.2+1.9)", and is possi- 
bly extended with an FWHM of (32.8±4.7)", which is a lower 
limit to its full width because of insufficient background signal. 
Indeed. lSciineider et all d2006l) find that the W75N cloud extends 
over several arcminutes. 



4. Discussion and Conclusions 

To estimate the H2O and 13 CO abundances from our data, 
we have run spherical radiative transfer models following 
Marseille et al. (2008). First, the d ust continuum emission was 
modeled with the MC3D program (Wolf & Henning 2000) with 
the source size and luminosity kept fixed at the values in iQ] The 
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Fig. 3. Plots of observed intensity versus spatial offset with 
Gaussian models superposed. 



continuum data are consistent with a power-law density profile 
n - rio(r/roy a with the index a- 1 .5 as expected for evolved pro- 
tostellar envelopes dVan der Tak et al. 2000). Derived tempera- 
tures range from 117K at the adopted inner radius of 0.01 pc 
to 23 K at the outer radius of 0.3 pc; densities drop from 
3xl0 7 cirr 3 to 2xl0 5 cirT 3 . This temperature and density pro- 
file was adopted for the line radiative transfer with the RATRAN 
program (Hogerheiide & van der Tak 2000). The abundance of 
H2O was varied between 10" ° and 10~ 7 and the 13 CO abun- 
dance between 10~ 7 and 10~ 6 , both independent of radius. 

The red lines in Figure [4] show the results of our best-fit 
model. The fit to the line profiles at the central few positions 
is good if two components are added to the dense core model: 
one representing the outflow and one for the foreground cloud. 
This simple addition procedure is valid because the foreground 
cloud is transparent and the outflow is seen at an angle, so that 
neither component blocks the view of the dense core. The mod- 
eled optical depths are 2.8 for the H2O line and 0.77 for 13 CO, 
and the derived abundances are uncertain to «35% which is the 
quadratic sum of 30% from the line intensity and 20% from the 
core mass estimate. The model underproduces the observed line 
strengths at large position offsets, which suggests that the den- 
sity profile flattens out at large radii, as also indicated by the 
observed continuum brightness at the offset positions. 

Deriving the 13 CO and H2O abundances in the outflow and 
the foreground using radiative transfer models is not possible be- 
cause the masses and H2 column densities of these components 
are unknown. Instead we have used RADEX dVan der Tak et al.1 
2007) to estimate the ir 13 CO and H2O c olumn densities, and use 
a 12 C/ 13 Cratioof 60 (M ilam et al.1l2005l) and a CO abundance of 
2xl0~ 4 dLacv et al.ll 1994 to convert 7Y( 13 CO) to N(U 2 ) and es- 
timate jc(H20) as 7Y(H20)/7Y(H2). For the foreground cloud, we 
adopt rkin= 10 K and n(H2)= 10 4 crrT 3 , and for the outflow, we 
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Fig. 4. Data from Fig. [2] with radiative transfer models super- 
posed. 

Table 1. Abundances and column densities of 13 CO and H2O. 



Component N( n CO) a JV(p-H 2 0) 2V(H 2 ) *(p-H 2 0) 

10 16 cm- 2 10 12 cm- 2 10 21 cm" 2 



Dense core 7.8xl(r 7 

Outflow 5 10,000 

Foreground 0.7 4 



1.6X10" 10 
15 7xl0- 7 

2.1 4x10-" 
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": Abundance for the dense core 



adopt 7\i n = 200K and «(H2)= 3xl0 4 cm- 3 , which assumptions 
introduce a factor of ~2 uncertainty each. To estimate Af( 13 CO) 
for the foregroun d cloud, we use the 7=1-0 observations by 
iJakob et al.l (I2007I) . 

Table Q] summarizes our derived column densities and abun- 
dances of 13 CO and H2O in the various physical components 
of the DR21 region. Our H2O abundance in the dense core is 
~100x lower than previous determinations (§Q3 but should be 
regarded as a lower limit. At the low temperatures and high den- 
sities in the core, most H2O is likely frozen on grains, and the 
observed line may arise in a small region with a high H2O abun- 
dance. The derived I3 CO abundance for the core is ~4x lower 
than expected for the above values of the CO isotopic ratios and 
abundance, which suggests that even some CO is frozen out in 
the outer parts of the core. The density of the foreground cloud 
is too low for significant freeze-out, but with A y a 1.2 mag, pho- 
todissociation is rapid for H2O but not for 13 CO. 

The high H2O abundance for the outflow is likely related to 
its temperature of ~200K, which is high enough to have H2O 
released from t he dust grains by th ermal evaporation, or possi- 
bly by shocks (Mel nick et al.1120081) . Further enhancement may 
be expected in even warmer gas ( ; 250 K) when neutral-neutral 
reactions drive most gas-phase oxygen into H2O, but such gas 
is not probed by our data. Future HIFI observations of high- 
excitation H2O lines towards protostars of all masses will how- 
ever very likely reveal this effect. 
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